Detection of two-sided alternatives in a Brownian Motion mode 
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Introduction and motivation 

The need for statistical surveillance has been noted in many different areas, including quality 
control (see for example [2]), epidemiology (see for example [13] ) , medicine (see for example [1]), 
machinery monitoring, seismology, finance (see for example [I]) etc. In this work, we address the 
problem of the detection of two-sided alternatives in a Brownian motion model. This model is the 
continuous time equivalent to the discrete time Gaussian observation model. For stochastic systems 
with linear dynamics and linear observations that are driven by Gaussian noise, the Kalman-Bucy 
innovation process is known to be a sequence of independent Gaussian random variables. Such models 
Q\ \ can be used to study systems subject to system component failures and other systems involving small 
non-linearities ( |161 I12j). Fault detection in a navigation system, where an abrupt change in the 
model parameters corresponds to an abrupt change in the mean of the Kalman filter innovations is an 
instance of such a situation ( JTOj, [2]). The sign of the change depends on the signs of the gyro errors 
(|15j). Another instance of such a model can be seen in sensor failure detection for the monitoring of 
traffic incidents on freeways. Each sensor is placed in different locations on the freeway and records 
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the mean velocity and density of cars. An abrupt and systematic change in these recordings would 
trigger an abrupt change in the Kalman filter innovations in either direction depending on whether 
the sensor is consistently overestimating or underestimating (see p3]). Identification and removal of 
the faulty sensor becomes essential. The continuous version of the Kalman filter innovations in all of 
the above linear Gaussian models is seen to be a Brownian motion ([8]). Other applications includes 
the detection of a rhythm jump of the heartbeat during an ECG (see [3]) and in the detection of a 
positive or negative drift in the log of stock price dynamics. 

This paper is concerned with the quickest detection of two-sided alternatives in the drift of 
a Brownian motion. In particular, we find the best 2-CUSUM stopping rule with respect to an 
extended Lorden criterion. Although, the mathematical formulation is done in the context of the one- 
dimensional case, extension to the vector case that corresponds to the Kalman innovations in linear 
systems described above is straightforward (see [7]). 



lr This research was supported in part by the U.S. National Science Foundation under Grant No. CCR-02-05214. 



Mathematical formulation and main results 



We sequentially observe a process {it} with the following dynamics: 
dwt t < T 



\i\dt + dw t 
or 

t -fi 2 dt + dw t 



t > T 



where r, the time of change, is assumed to be deterministic but unknown; wt is a standard Brownian 
motion process; fa, the possible drifts to which the process can change, are assumed to be known, 
but the specific drift to which the process is changing is unknown. Both [i\ and fi 2 are assumed to be 
positive. 

The probability triplet consists of (C[0, oo], U$>o.Ft)> where Tt = er{£ s ,0 < s < t} and the families 
of probability measures {P*}, r G [0, oo), whenever the change is fa, i = 1,2, and Voo, the Wiener 
measure. 

Our goal is to detect a change by means of a stopping rule T adapted to the filtration T%. As a 
performance measure for this stopping rule we propose an extended Lorden criterion (see [5]) 

(1) Jl{T) = maxsup essup E l T (T — t) + \!F t . 

This gives rise to the following min-max constrained optimization problem: 



(2) 



inf J L {T) 
subject to £"00 [T] > 7, 



where the constraint specifies the minimum allowable mean time between false alarms. 

In this paper we seek the best 2-CUSUM stopping rule in the sense described in (|2|). The 2-CUSUM 

rules have been proposed and used extensively due not only to the simplicity in the calculation of 

their first moment (see [9]), but also to their asymptotically optimal character (see [5], [TT]). 

We begin by defining the CUSUM statistics and stopping rules of interest. 

Definition Let v\ > and v 2 > 0. Define 
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it - 5fJ.1t; mf = inf s < t u+; yf 
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<t - 2fi2t; m~ t 



inf 



s<t u s ; Vt 
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3. Ti(vi) = inf{t > 0; yf > ui}, and 

4. T 2 {v 2 ) =inf{t > 0;yt > ^2}. 

The 2-CUSUM stopping rules are then of the form T(v\, v 2 ) = T\{v{) A T 2 (v 2 ). 

We also define the following stopping rules, the use of which will become apparent later. 

Definition For a > and b > 0, we define 

1. U+(a) = ia£{t > 0; uf > a}, 

2. I/- (6) = inf{t > 0; -uj < -b}, and 

3. n(a, b) = P (U + (a) < U~(b)) . 



For any 2-CUSUM stopping rule T we have (see [5]) J L (T) = max{E% [T],E$ [T]}. 
We now classify 2-CUSUM rules according to the class Q = {T{ui,v 2 );vi = z^} of harmonic mean 
rules and the classes C\ = {T{v\, v 2 ) \ v\ > u 2 > 0} and C 2 = {T{y\, v 2 ) I v 2 > v\ > 0} of non-harmonic 
mean rules. For simplicity of display and notation we finally define the constants m = min{i/i, i^}, 
M = max{i^i,i/ 2 } and the functions C m (x,y) = /m( i™+j m (y) ' X x(y) = y f x ( y ) +x > fy( x ) = fx(y) = 

e y —yx—i summarize the main results. 

!r 

Theorem Let T{y\, 1*2) = Ti(yi) A T 2 {y 2 ) be any 2-CUSUM stopping rule and denote T(u\, u 2 ) by T. 
Then, the following is true under any of the measures Poo, Pq 1 and P 2 : 



1. for all T £ C\, m = V2, M = u\, we have 



E[T] = E[T 2 (m)} 



/ 1 ^ (M—m)n 

, — lim II ( — , m 

E[T 1 (m)}+ E[T 2 {m)]n^oo \ n ' 



E[T 2 (m)} 



and 



2. for all T E C 2 , m = v\, M = u 2 , we have 
E[T) = S[Ti(m) 
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Corollary Let T(Vi,f 2 ) = T\{y\) f\T 2 [y 2 ) be any 2-CUSUM stopping rule and denote T(vi,v 2 ) by 
T. Then, for all T G Ci, m = i/ 2 , M = v\ and 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



Eoc[T] < 2/ m (/i 2 ) 

Poo PI > 2/ m 0ia) 



1 _ C m (n 2 ,m) c _ Am (-^)(M-m) 
/m(M2) 

1 _ C m (n2,Hl) c -\ m (u 2 )(M-m) 



E l [T] < 2/ m (^ 2 + 2/ii 

PfJPI > 2/ m (/x 2 + 2 m 

PqPI < 2/ m (-/x 2 ) 

PqPI > 2/ m (-/x 2 ) 



j _ Cm(^2 + 2^1, -/ij) / -A ro (^!)(M-m) 

/m(^2 + 2/ij) 
1 _ C m{^2 + 2/ii, ^/fl) e -A m (jU 2 +2jHl)(M-m) 
/to(M2 + 2 Ml) 
1 _ C m (-/i 2 ,/Xl + 2/X 2 ) c _A m (-( Atl+ 2 M2 ))fM-m) 



/m(— ^2) 

C m (-/i 2 ,/xi + 2^ 2 ) 



, and 



-A m (-^2)(Af-m) 



/m(— ^2) 

Similar results hold for T € C 2 . For more details please refer to [6]. 

Theorem The best T* 2-CUSUM stopping rule exists and is unique and we distinguish the following 

cases 

1. If iii < fi 2 then T* G C 2 . 

2. If fi 2 < Ml then X* 6 Ci- 

3. If mi = M2 then T* G Q. 

We refer the reader to [6] for a detailed proof of all of the above results and other interesting corollaries. 



References 



[1] E. Andersson. Monitoring cyclical processes: A non-parametric approach. Journal of Applied 
Statistics, 29:973-990, 2002. 

[2] M. Basseville and I. Nikiforov. Detection of Abrupt Changes: Theory and Applications. Prentice 
Hall, Englewood Cliffs, NJ, 1993. 

[3] P.C. Doerschuk, R. R. Tenney, and A. S. Willsky. Estimation-based Approaches to Rhythm 
Analysis in Electrocardiograms, Lecture Notes in Control and Information Science III, pages 297- 
313. Springer- Verlag, New York, 1986. 

[4] M. Frisen. Evaluations of methods for statistical surveillance. Statistics in Medicine, 11:1489- 
1502, 1992. 

[5] O. Hadjiliadis and G. V. Moustakides. Optimal and asymptotically optimal CUSUM rules for 
change point detection in the Brownian motion model with multiple alternatives. Theory of 
Probability and Its Applications (Teoriya Veroyatnostei i ee Primeneniya), 50(1):131-144, 2006. 

[6] O. Hadjiliadis and H.V. Poor. The best 2-CUSUM stopping rules for quickest detection of two- 
sided alternatives in a Brownian motion model. 2007. submitted. 

[7] G. V. Moustakides. Optimality of the CUSUM procedure in continuous time. Annals of Statistics, 
32(1):302-315, 2004. 

[8] B. Oksendal. Stochastic differential equations. Springer- Verlag, Berlin, 1985. 

[9] D. Siegmund. Sequential Analysis. Springer- Verlag, New York, 1st edition, 1985. 

[10] A. Tartakovsky. Minimal time detection algorithms and applications to flight systems. December 
1993. Technical report, Flight Systems Research Center, University of California, Los Angeles. 

[11] A. G. Tartakovsky. Asymptotically minimax multi-alternative sequential rule for disorder de- 
tection. Statistics and Control of Random Processes: Proceedings of the Steklov Institute of 
Mathematics, 202(4) :229-236, 1994. American Mathematical Society, Providence, RI. 

[12] R. R. Tenney, R.S. Hebbert, and N. R. Sandell. A tracking filter for maneuvering sources. IEEE 
Transactions on Automatic Control, 22(2):246-251, April 1977. 

[13] J.A.C. Weatherall and J.C. Haskey. Surveillance of malformations. British Medical Bulletin, 
32:39-44, 1976. 

[14] A.S. Willsky, E. Chow, S. B. Gershwin, C. S. Greene, P.K. Houpt, and A. L. Kurkjian. Dy- 
namic model-based techniques for the detection of incidents on freeways. IEEE Transactions on 
Automatic Control, 25(3):347-359, 1980. 

[15] A.S. Willsky, J.J. Deyst, and B.S. Crawford. Two self-test methods applied to an inertial system 
problem. Journal of Spacecraft, 12(7):434-437, 1975. 

[16] A.S. Willsky and H. Jones. A generalized likelihood ratio approach to the detection and estimation 
of jumps in linear systems. IEEE Transactions on Automatic Control, 21(1):108-112, February 
1976. 



